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Abstract 

We study the problem of assigning transmission ranges to radio stations placed arbitrarily in 
a d-dimensional (d-D) Euclidean space in order to achieve a strongly connected communication 
network with minimum total power consumption. The power required for transmitting in range r 
is proportional to r“, where a is typically between 1 and 6, depending on various environmental 
factors. While this problem can be solved optimally in ID, in higher dimensions it is known to be 
AP-hard for any a > 1. 

For the ID version of the problem, i.e., radio stations located on a line and a > 1, we propose 
an optimal 0(n^)-time algorithm. This improves the running time of the best known algorithm by 
a factor of n. Moreover, we show a polynomial-time algorithm for finding the minimum cost range 
assignment in ID whose induced communication graph is a t-spanner, for any t > 1. 

In higher dimensions, finding the optimal range assignment is AP-hard; however, it can be 
approximated within a constant factor. The best known approximation ratio is for the case a = 1, 
where the approximation ratio is 1.5. We show a new approximation algorithm with improved 
approximation ratio of 1.5 — e, where e > 0 is a small constant. 


1 Introduction 

A wireless ad-hoc network is a self-organized decentralized network that consists of independent radio 
transceivers (transmitter/receiver) and does not rely on any existing infrastructure. The network nodes 
(stations) communicate over radio channels. Each node broadcasts a signal over a fixed range and 
any node within this transmission range receives the signal. Communication with nodes outside the 
transmission range is done using multi-hops, i.e., intermediate nodes pass the message forward and 
form a communication path from the source node to the desired target node. The twenty-first century 
witnesses widespread deployment of wireless networks for professional and private applications. The field 
of wireless communication continues to experience unprecedented market growth. For a comprehensive 
survey of this field see [7j. 

Let 5" be a set of points in the d-dimensional Euclidean space representing radio stations. A range 
assignment for S' is a function p : S —>■ R''" that assigns each point a transmission range (radius). 
The cost of a range assignment, representing the power consumption of the network, is defined as 
cost[p) = F)r some real constant a > 1, where a varies between 1 and values higher than 

6, depending on different environmental factors [7]. 

A range assignment p induces a directed communication graph Gp = {S,Ep), where Ep = {(m,u) : 
p{u) > |uu|} and |mu| denotes the Euclidean distance between u and v. A range assignment p is valid if 
the induced (communication) graph Gp is strongly connected. For ease of presentation, throughout the 
paper we refer to the terms ‘assigning a range \uv\ to a point u G S’ and ‘adding a directed edge (m,u)’ 
as equivalent. 

We consider the d-D Minimum Cost Range Assignment (MinRange) problem, that takes as 
input a set S' of n points in R'^, and whose objective is finding a valid range assignment for S of 
minimum cost. This problem has been considered extensively in various settings, for different values of 
d and a, with additional requirements and modifications. Some of these works are mentioned in this 
section. 
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Ini, Kirousis et al. considered the ID MinRange problem (the radio stations are placed arbitrarily 
on a line) and showed an 0(n'^)-time algorithm which computes an optimal solution for the problem. 
Later, Das et al. [5] improved the running time to 0(n^). Here, we propose an 0(n^)-time exact 
algorithm, this improves the running time of the best known algorithm by a factor of n without increasing 
the space complexity. The novelty of our method lies in separating the range assignment into two, left 
and right, range assignments (elaborated in Section . This counter intuitive approach allows us to 
achieve the aforementioned result and, moreover, to compute an optimal range assignment in ID with 
the additional requirement that the induced graph is a t-spanner, for a given t > 1. 

A directed graph G = {S, E) is a t-spanner for a set S, if for every two points u,v £ S there exists a 
path in G from m to n of length at most t|un|. The importance of avoiding flooding the network when 
routing, was one of the reasons that led researchers to consider the combination of range assignment and 
t-spanners, e.g., [Dimn], as well as the combination of range assignment and hop-spanners, e.g., BB- 
While bounded-hop spanners bound the number of intermediate nodes forwarding a message, t-spanners 
bound the relative distance a message is forward. For the ID bounded-hop range assignment problem, 
dementi et al. [1] showed a 2-approximation algorithm whose running time is 0(hn^). To the best of 
our knowledge, we are the first to show an algorithm that computes an optimal solution for the range 
assignment with the additional requirement that the induced graph is a t-spanner. 

While the ID version of the MinRange problem can be solved optimally, for any d > 2 and a > 1, 
it has been proven to be A^P-hard (in for d > 3 and 1 < a < 2 and later in [3] for d > 2 and 
a > 1). However, some versions can be approximated within a constant factor. For a = 2 and any d > 2 
Kirousis et al. [5] gave a 2-approximation algorithm based on the minimum spanning tree (although they 
addressed the case of d G {2,3} their result holds for any d > 2). The best known approximation ratio 
is for the case a = 1, where the approximation ratio is 1.5 [2]. We show a new approximation algorithm 
for this cas^with improved approximation ratio of 1.5 — e, for a suitable constant e > 0. We do not 
focus on increasing e but rather on showing that there exists an approximation ratio for this problem 
that is strictly less than 1.5. This is in contrast to classic problems, such as metric TSP and strongly 
connected sub-graph problems, for which the 1.5 ratio bound has not yet been breached. 

2 Minimum Cost Range Assignment in ID 

In the ID version of the MinRange problem, the input set S = {vi,...,Vn] consists of points located 
on a line. For simplicity, we assume that the line is horizontal and for every i < j, Vi is to the left of Vj. 
Given two indices 1 < i < j < n, we denote by Sij the subset {vi, ..., Vj} C S. 

We present two polynomial-time algorithms for finding optimal range assignments, the first, in Sec¬ 
tion (2?^ for the basic ID MinRange problem, and the second, in Section [2?^ subject to the additional 
requirement that the induced graph is a t-spanner (the ID MinRangeSpanner problem). Our new 
approach for solving these problems requires introducing a variant of the range assignment. Instead of 
assigning each point in 5" a radius, we assign each point two directional ranges, left range assignment, 
p’’ : S ^ R'*', and right range assignment, p'^ : S ^ R+. A pair of assignments (p^p’") is called a 
left-right assignment. Assigning a point v G S a, left range p*(n) and a right range p^{v) implies that 
in the induced graph, Gpir-, v can reach every point to its left up to distance p\v) and every point to 
its right up to distance p''{v). That is, Gpir-, contains the directed edge {vi,Vj) iff one of the following 
holds: (i) i < j and \viVj\ < p^{vi), or (ii) j < i and \viVj\ < p^{vi). The cost of an assignment (p^,p’'), 
is defined as cost{p'‘, p'') = 

Our algorithms find a left-right assignment of minimum cost that can be converted into a range 
assignment p with the same cost by assigning each point v G S a range p{v) = max{p^(w), p’'('(;)}. Note 
that any valid range assignment for S can be converted to a a left-right assignment with the same cost, 
by assigning every point v G S, p^{y) = p’'(n) = p{v). To be more precise, either p^{v) or p’'(n) should be 
reduced to |um|, where u is the farthest point in the directional range (for Lemmato hold). Therefore, 
a minimum cost left-right assignment, implies a minimum cost range assignment. 

In addition to the cost function, we define cost'{p^, p'') = refine the 

term of optimal solution to include only solutions that minimize cost'{p^, p*") among all solutions, (p^, p*"), 
with minimum cost{p^, p''). 

^Values of a smaller than 2 correspond to areas, such as, corridors and large open indoor areas [7]. 
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2.1 An Optimal Algorithm for the ID MinRange Problem 


Das et al. [5] state three basic lemmas regarding properties of an optimal range assignment. The following 
three lemmas are adjusted versions of these lemmas for a left-right assignment. 

Lemma 1. In an optimal solution {p^p^) for every Vi G S, either p\vi) = 0 or p^{vi) = \viVj\ and 
similarly, either p^{vi) = 0 or p'"{vi) = \viVk\ for some j < i < k. 


Lemma 2. Given three indices l<i<j<k<n, consider an optimal solution for Si^k, denoted hy 
{p^,p'^), subject to the condition that p^{vj) > \viVj\ and p^{vj) > \vjVk\, then, 

• for all m = i, ...,j - I, p^{vra) = |rmWm+i| and p\vm) = 0; and 

• for all TO = j + 1, ...,k, p’-{Vm) = \VjnVm-l\ and p^{Vm) = 0. 


Lemma 3. In an optimal solution {p\p^), p^{vi) = 0 and p’'(ui) = |wir 2 |. 

Lemmaallows us to simplify the notation p^{vi) = \viVj\ for x € {l,r} and I < i,j < n, and write 
p^{i) = j for short. We solve the MinRange problem using dynamic programming. Given 1 < i < n, we 
denote by OPT{i) the cost of an optimal solution for the sub-problem defined by the input Si^m subject 
to the condition that p’'(z) = i -I- 1. Note that the cost of an optimal solution for the whole problem is 
OPT(l). 

In Section |2.1.l| we present an algorithm with 0{n^) running time and 0{n^) space (the same time 
and space as in [^j. Then, in Section 2.1.2 we reduce the running time to 0{n^). 


2.1.1 A Cubic-Time Algorithm 

Algorithm IDMinRA (Algorithm applies dynamic programming to compute the values OPT(i) for 
every 1 < z < n and store them in a table, T. Finally, it outputs the value T[l]. In our computation 
we use a 2-dimensional matrix. Sum, storing for every 1 < z < j < n the sum J2m=i While 

Algorithm 1 IDMinRA(S') 
for z = n — 1 downto 1 do 
for j = n downto z -|- 1 do 


Sum[z, j] ^ 

for z = rz — 1 downto 1 do 


kn-W«|“ ,z = n-l 

Sum[z-I-l,n]-I-luiUi+il" ,j = n 
Sum[z, j -I- 1] — \vjVj+i\°' ,otherwise 


T[z]^ 


2|z;iUi+i|“ 

min {Sum[z, fc' — 1] -|- T[fc' — 1] — \vk'-iVk' 

i<.k<n 
k<k' <n 

max{\viVk\°', lufeUfc'T}} 


, z = n — 1 
, otherwise 


return T[l] 


the table T maintains only the costs of the solutions, the optimal assignment can be easily retrieved by 
backtracking the cells leaded to the optimal cost and assigning the associated ranges (described in the 
proof of Lemma . 


Correctness. We prove that for every 1 < z < n, the value assigned to cell T[z] by the algorithm equals 
OPT(i). Trivially, OPT(n— 1) indeed equals 2\viVi+i\°‘. Assume, during the z-th iteration it holds that 
T[z'] = OPT{i') for every i < i' < n, the correctness of the computation done during the z-th iteration 
is given in Lemma 

Lemma 4. Given an index i with 1 < i < n — 1, 

( k'-2 

OPT{i) = min I ^ \vmVm-i-i\°‘+ OPT{k'- 1) - \vk>-iVk’\°'+ T[nax{\viVk\°',\vkVk’\°‘} 

k<k^<n \ rn—i 
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Proof. Let Xi denote the right side of the equation, we prove OPT{i) = Xi. 

OPTii) < Xi'. We show that all costs that appear as min function arguments in Xi correspond to valid 
assignments and thus infer, by the optimality of OPT{i), that the above inequality holds. Consider an 
argument with parameters k and k'. We associate it with an assignment (p^p^) defined as follows (see 
Fig. [^a)). For m>k' — l the assignment is inductively defined by OPT{k' — 1). For every i < m < k, 
p^{m) = m and p'^{m) = m + 1, for every k < m < k', p\m) = m — 1 and p'^{m) = m {vk'-i is 
reassigned) and for k, p^{k) = i, p'^{k) = k'. By the validity of OPT{k' — 1), every two points among 
Sk'-i,n are (strongly) connected. Our assignment for Si^k'-i guarantees the connectivity between every 
two points in Si^k>-i, and thus between every two points in Si^n- 

OPT{i) > Xp. Consider an optimal solution (p^p’") for the points Si^n subject to the condition that 
p''{i) = i + 1. Let Vk be a point to the right of Vt with p*(fc) = i and let p''(fc) = k'. Note that since Vi 
is the leftmost point and the induced graph is strongly connected, such a point necessarily exists. 

Next we show that there is no edge directed either left or right connecting two points on different 
sides of Vk' in Gpir, except for possibly an edge {vj,Vk>-i) with j > k'. Assume towards contradiction 
that the former does not hold, i.e., there exists i < t < k', with p’’(t) > k'] then, reassigning p^{k) = 
max{t, k} maintains the connectivity, and reduces the value of cost' without increasing the value of 
cost in contradiction to the optimality of the solution. Now, let Vj be a point to the right of Vk' with 
p\j) = / G [b 'we show that j' > k' — 1. Consider a point j' < t < k', as we have shown, p’’(t) < k'. 
By symmetric arguments we have p\t) > f (see Fig. [^b)). Namely, there is no edge going out of 
the interval {vj',Vk'). Thus, connectivity can be achieved only if this interval is empty of vertices, i.e., 
either f = k' — 1 or f = k' (note that k' — \ > i). 

The above observation allows us to divide the problem into two independent subproblems, one for 
the points Si^k'-i subject to the constraints p\k) = i and p’'(fc) = fc', and the other for the points 
Sk'-i,n subject to the artificial constraint p''{k' — 1) = k' that guarantees the existence of a path 
from fc' — 1 to A:', due to the solution of the first subproblem, but should not be paid for. Regarding 
the first subproblem, by Lemma in an optimal assignment, for every i < m < k, p^{m) = m and 
p''im) = m + 1, and for every k < m < k' — 1, p^{m) = m — 1 and p''{m) = m. Thus, its cost is 
'^Tn=i + max{|ufcUi|“, \vkVk'\°‘'\. The cost of an optimal solution to the second subproblem 

is OPTik' — 1) — Hence, the cost of an optimal solution to the whole problem is the sum 

of the above costs and the lemma follows. 


□ 



o- o 



o o o 


o 

v„ 


(b) 


Figure 1: (a) An illustration of the assignment associated with OPT{i) with respect to given parameters 
k and k'. In gray are range assignments associated with OPT{k' — 1). (b) An illustration of the proof of 
Lemma 1^ In dashed arrows, the impossible ranges of Vt and in gray, the alternative assignment of lower 
cost'. 


Complexity. Obviously, Algorithm IDMinRA requires 0{'nf ) space. Regarding the running time, 
0{n) iterations are performed during the algorithm, each iteration takes 0{nf ) time Therefore, the total 
running time is 0{n^) and Lemma [^follows. 

Lemma 5. Algorithm IDMinRA runs in 0{n^) time using 0{n^) space. 

2.1.2 A Quadratic-Time Algorithm 

In this section we consider Algorithm IDMinRA from previous section and reduce its running time to 
O(n^). Consider the equality stated in Lemma Observe that given hxed values i and /c', the value 
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k that minimizes the argument of the min function with respect to i and k' is simply the value k that 
minimizes max{|z;iUfe|“, \vkVk'\°‘}- This value is simply the closest point to the midpoint of the segment 
ViVk', denoted by c(i,k'). Thus, 


OPT(i) = min 

2+l<A;'<n 


|u„Um+i|“ + OPT{k' - 1) - \vk'-iVk'\°‘ 

m—i 

+ max{|wjWc(i,fe-)|“, \vc(i,k')Vk'r} 


Consider Algorithm IDMinRA after applying the above modification in the computation of T[i\. Since 
there are only 0{ri) sub-problems to compute, each in 0(n) time, the running time reduces to 0{n?) 
and the following theorem follows. 

Theorem 6. The ID MinRange problem can be solved in 0(n^) time using O(n^) space. 


2.2 An Optimal Algorithm for the ID MinRangeSpanner Problem 

Given a set S = {wi,u„} of points in ID and a value t > 1, the ID MinRangeSpanner problem aims 
to find a minimum cost range assignment for S, subject to the requirement that the induced graph is a 
t-spanner. We present a polynomial-time algorithm which solves this problem optimally and follows the 
same guidelines as Algorithm IDMinRA. 

We begin with providing the key notions required for understanding the correctness of the algorithm, 
followed by its description. Due to space limitation, we do not supply a formal proof. The first and 
most signihcant observation, is that the problem can still be divided into two subproblems in the same 
way as in Algorithm IDMinRA, by similar arguments to those of Lemma In Lemma we show 
that any assignment that does not satisfy the conditions required for the division can be adjusted to a 
new assignment with a lower value of cost' that preserves connectivity. The new assignment, however, 
preserves also the lengths of the shortest paths, which make the argument legitimate for this problem as 
well. 

The two problems (MinRange and MinRangeSpanner ) differ when it comes to solving each of 
the above subproblems. Consider the left subproblem, i.e., of the form described in Lemma The 
optimal assignment for it is no longer necessarily the one stated in the lemma, since it does not ensure 
the existence of t-spanning paths. Therefore, our algorithm divides problems of this form into smaller 
subproblems handled recursively (see Fig. right). Dealing with such subproblems, requires defining 
new parameters: a rightmost input point Vj, and the length of the shortest paths connecting Vi to Vj, 
Vj to Vi and Vi to Vi+i not involving points in Sij except for the endpoints, denoted by , 6 , and < 5 * , 
respectively. Regarding the computation of a subproblem, since points may be covered now by vertices 
outside the subproblem domain, we allow Vk to have either a right or a left range equals 0 (in the terms 
of Algorithm IDMinRA, either k = i oi k = k'). 

Another key observation is that any directed graph G over S' is a t-spanner for S iff for every 1 < i < n 
there exists a t-spanning path from vt to Vi+i and from to Vi. Moreover, given that G is strongly 
connected implies that the addition of an edge between consecutive points does not effect the length of 
the shortest path between any other pair of consecutive points. Therefore, for subproblems with j = i + 1 
we assign p’’(i) = i -I-1 (resp. p^{i -I- 1) = i) iff /\i, i + \\ > t (resp. ^/|z, i -|- 1| > t) and thus ensuring 
that the induced graph is a t-spanner. 

Our algorithm may consider solutions in which an assignment to a node is charged more than once 
in the total cost; however, for every such solution, there exists an equivalent one in which the charging 
is done properly and is preferred by the algorithm due its lower cost. 

We denote by OPT{i,j, ^, S , S') the cost of an optimal solution to the sub-problem defined by the 
input Sij subject to the parameters (5 , and < 5 * representing the lengths of the shortest external paths 
as defined earlier in this section. Let Ai j = {2\viVi\ -I- \viVj\, 2\vjVr\ + \viVj\ '■ I < i < j < r}. We compute 
OPT{iG,~^,S') for every 1 < i < j < n, and the corresponding € Aij iteratively, while in 

stage X all subproblems with j — i = x are solved. The computation is derived from the equalities below. 
For simplicity of presentation, we overload notation and write \i,j\ to mean \viVj\. In addition, we write 
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0 in place of (5* where (5* = "?. 

OPT{i, z + 1,(5®) = + V, where 

/ I*, * + 1| (*assj3riOT5 p®'(z) = z + 1*) , "^/jz, z + 1| > f 


, otherwise 


and V is defined symmetrically. For j > 1, we have 

opr(z,j,t,t,d®) = 

|z, z + 1|“ + OPT{i, z + 1, |z, z + 1|, |z + z, j| + d , 0) 

+ OPT{i + 1, j, oo, (5 — |z + 1,z|, 0) 

|z,z + ir+OPr(zz + l, d®, |f,z + l|, 0) 

+ OPT{i + 1, j, |z,z + 1| + "?, d — |z,z + 1|, 0) 

max{|z, fc|“, \k, k'\°‘} 

+ OPT{i, z + 1, 5\ |z + 1, A:| + |fc, z|, 0) 

+ OPT{i + 1, k, oo, |A;,z + l|, oo) 

+ OPT{k, k' — 1, |fc, k' — 1|, oo, \k, k + 1|) 

+ OPT{k' - 1, j, \k' - 1, z| + 7, Y - |z, k' - 1|, \k' -l,k\ + \k, k'\). 


min < 

i^k<j 
k<k' <j 


,i = k 
= k' 

,i ^ k k' 


We permit either i = k and then fc® = z + 1 or fc = fc' and then A:' = z + 1 but not both. 


I I 
9 —►tp o o 

Vi W+1 


V, 


i = k and k' = i + 1 


I I 

—9 o o 

W W+i 


k = k' = i + 1 




\Ar^ Xl/"''" ^ 

(f 90 09 009 009 

Vk ~ Vy Vj 

i ^ k ^ k' 


Vi 


Figure 2: An illustration of the algorithm for the MinRangeSpanner problem. The ranges are illus¬ 
trated in black arrows and the division to subproblems in sashed lines. 


Complexity. Let A be the set of all distinct distances in S, then for every Vi,Vj G S, |Aij| = |A| = 
0{n). We fill a table with 0(n^|Ap) cells, each cell is computed in 0{'n?) time, thus, the total running 
time is 0(n^|Ap). As we have focused on presenting a simple and intuitive solution, rather than reducing 
the running time, a more careful analysis achieves a better bound on the time complexity. For example, 
the relevant domain of ~i, 6 , and d® can be estimated more precisely with respect to t. Moreover, 
Observation allows reducing the running time by a factor of n. This is done by decreasing the number 
of relevant combinations of z and k' that have to be checked by the algorithm, for fixed indices j and k 
with i < k < k' < j, to 0(n), using similar arguments to those in Lemma|^ (see Fig. [^. 

Observation 7. Consider an optimal assignment (p*,p®’) and a point Vk G S. Let p^{k) = z and let fc® 
denote the minimal index with k < k^ and IvkVyl > \vkVi\, then /c®+^ < p^{k) < fc®. 



O /O o 

Vf^i '' 


o o 

Vn 


Figure 3: An illustration of Observation!^ Every pair of symmetric arcs indicates equal distances from 
Vk- The marked domains indicates the legal values of p^{k) for different values of z. 
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3 The MinRange Problem in Higher Dimensions 

In this section we focus on the MinRange problem for dimension d> 2 and a = 1. As all the versions 
of the problem for d > 2 and a > 1, it is known to be A^P-hard. Currently, the algorithm achieving the 
best approximation ratio for a = 1 and ant d > 2 is the Hub algorithm with a ratio of 1.5. This algorithm 
was proposed by G. Calinescu, P.J.Wan, and F. Zaragoza for the general metric case, and analyzed by 
Ambiihl et al. in [5] for the restricted Euclidean case. We show a new approximation algorithm and 
bound its approximation ratio from above by 1.5 — e for e = 5/10®. Although in some cases our phrasing 
is restricted to the plane, all arguments hold for higher dimensions as well. 

In our algorithm we use two existing algorithms, the Hub algorithm and the algorithm for the ID 
MinRange problem introduced by Kirousis et al. [5], to which we refer as the ID RA algorithm. We 
observe that the later algorithm outputs an optimal solution for any ordered set V = {ui, ...,u„} with 
distance function h that satisfies the following line alike condition: for every I < i < j < fc < / < n, it 
holds that h{vi,vi) > h{vj,Vk)- We use this algorithm for subsets of the input set that roughly lie on a 
line. 

3.1 Our Approach 

Presenting our approach requires acquaintance with the Hub algorithm. The Hub algorithm finds the 
minimum enclosing disk C oi S centered at point hub G S. Then, the algorithm sets p{hub) = rmin, 
where is C’s radius. Finally, it directs the edges of MST{S) towards the hub and for each directed 
edge {v,u) sets p{v) = \vu\. The cost of this assignment is w{MST{S)) + rmin < w{MST{S)) + 
{w{MST{S)) + w{eM))/‘2, where Cm is the longest edge in MST{S) and the weight function w is defined 
with respect to Euclidean lengths. 

To guide the reader, we give an intuition and a rough sketch of our algorithm. We characterize the 
instances where the Hub algorithm gives a better approximation than 1.5, and to generalize these cases 
we slightly modify it. Furthermore, we show an algorithm that prevails in the cases where the modified 
Hub algorithm fails to give an approximation ratio lower than 1.5. Before we elaborate more on the 
aforementioned characterization, another piece of terminology. Given a graph G over S and two points 
p,q G S, the stretch faetor from p to g in G is Scip, q)/\pq\, where Scip, q) denotes the Euclidean length 
of the shortest path between p and q in G. We use ^large when referring to values greater than fixed 
thresholds, some with respect to w{MST{S)), defined later. 

Consider MST{S) and its longest path Pm- If one of the following conditions holds, then the Hub 
algorithm or its modification results in a better constant approximation than 1.5: (Al) there exists a 
^large edge in MST{S); (A2) a ^large fraction of Pm consists of disjoint sub-paths connecting pairs of 
points with ^^large stretch factor, not dominated by one sub-path of at least half the fraction; or (A3) 
the weight w{MST[S)\Pm) is ^large. 

Otherwise, there are three possible cases: (Bl) the graph MST{S) is roughly a line; (B2) there are 
two points in Pm with '^large stretch factor, i.e., there is a ^large ‘hill’ in Pm, and then either MST{S) 
roughly consists of two ID paths; or (B3) the optimal solution uses edges connecting the two sides of 
the ‘hill’, covering ^large fraction of it. 

The last three cases are approximated using the following method. We consider every two edges 
connecting the two sides of the ‘hill’ as the edges in the optimal solution that separates the uncovered 
remains of the path to two independent sub-paths, i.e., not connected by an edge, (the points in both 
sub-paths may be connected to the middle covered area.) Note that such two edges exist. We direct the 
covered area to achieve a strongly connected sub-graph and solve each of the two sub-paths separately 
in two techniques. The first, using the \D RA algorithm with a distance function implied by the input, 
satisfying the line alike condition, and applying adjustments on the output, and the second, using the 
Hub algorithm. A (1.5 — e)-approximation is obtained for cases (Bl) and (B2), using the first technique, 
and for cases (Bl) and (B2), using the second technique. The algorithm computed several solutions, 
using the aforementioned methods, and returns the one of minimum cost. 

3.2 The Approximation Algorithm 

The algorithm uses the following three procedures that are defined precisely at the end of the algorithm’s 
description. 
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• The flatten procedure / - a method performing shortcuts between pairs of points on a given path P 
resulting in a path without two points of stretch factor greater than Cg. 

• The distance function hs - a, distance function defined for an ordered set PCS, satisfying the line 
alike condition. 

• The adjustment transformation g - a function adjusting an optimal range assignment for an ordered 
set PCS with distance function h, to a valid assignment for P. 

Let R be the forest obtained by omitting from MST{S) the edges of its longest path, PM. Given a 
point V G Pm, let T{v) denote the tree of R rooted at v. For every u G T{v) let r{u) denote the root of 
the tree in R containing u, namely, v. For a set of points V C Pm, let T{V) denote the union lj„gy T{v). 
For ease of presentation, we assume the path Pm has a left and a right endpoints, thus, the left and 
right relations over Pm are naturally defined. 


The main algorithm scheme: 

Compute four solutions and return the one of minimal cost. In case of multiple assignments to a point 
in a solution, the maximal among the ranges counts. 

Solution (i): apply the Hub algorithm. 

Solution (ii): apply a variant of the Hub algorithm - hnd a point c G Pm that minimizes the value 
Tc = max{|cpi|, |cpz|}, where pi and Pz are the endpoints of the path Pm- Assign c the range rc, direct 
Pm towards and bi-direct all edges in R. 

(* The rest of the algorithm handles cases (B1)-(B3) defined in Section 3.1 *) 

For every edge e G Pm do : 


Let Pel and Pg,- be the two paths of Pm\s, to the left and to the right of e, respectively. 
Apply the flatten procedure f on P^i and Pg,- to obtain the sub-paths 
Pi' = {pi,P2,-,Pm) and Pr' = {pjn+i,Pi-\-2,-■,Pz), respectively. 

(* Note R has been changed during the flatten procedure *) 

For every 4 points pppr ,pr' ,Pr with I < V < m < r' < r in the flattened sub-paths: 


In both solutions (iii) and (iv) direct the path P^ = {pi, ...pm,Pm+i, ■■■,Pr) towards pi and for 
each point pi with \ < i < z direct T(pi) towards pi and assign pi a range w{T(pi)). Perform 
the least cost option among the following two, either add the edge {pi,Pr), or add the two edges, 
one from u; to for ui G T{pi),Ur' G Tfpz-') of minimal length and the other from un to u^. for 
Ur G T{pi'),Ur G T(pr) of minimal length, (see illustration in Fig. [^. As for the two sub-paths 
Pi = {Pi,P 2 , ■■■,Pi) and Pr = {pr,Pr+i, ---tPz), assign them ranges as follows: 

Solution (iii): apply the Hub algorithm separately on each sub-path. 

Solution (iv): apply the IP RA algorithm separately on each sub-path with respect to the distance 
function hs and perform the transformation g on the assignment received. 



Figure 4: The two sub-paths Pi and Pr as defined in the algorithm. 


The flatten procedure /. Let Cs = 5/4. Given a path P = {ui,..,u„}, set Qp = {}. Let j > * be 
the maximal index such that Sp{vi,Vj) > Cs\viVj\. If such index does not exist, let j = i + 1. Else 
(j > i + 1), add the edge {vi,Vj) to P, remove the edge {vj-i,Vj) from P, move the sub-path 
{vi, ..,Vj-i) from P to the forest P, and update Qp = Qp U {{vi,Vj)}. Finally, repeat with the 
sub-path {vj,..,Vn) without initializing Qp. 

The definitions for hs and g are given with respect to the sub-paths Pi and Pr, the definitions for 
the sub-path Pr and Pp are symmetric. 
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The distance function / 15 . For every two points Pj,Pk with 1 < j < k < I we define, 

hs{Pj,Pk) = min |un|. 

ueT{p-,),i<j' <j 

vGT{pf^/),k<k' <171 


The adjustment transformation g. Given an assignment p' : Pi ^ R’*', we transform it into an 
assignment g{p') = p : Pp ^ R"*". First, we assign ranges a follows: 


p{Pj) = 


Cs ■ p{Pj) + Ck ■ T{pj), l<j<l, 
Cfe • T{pj), I <3 <m, 


where Ck = 1 + 8(1 + Cs) = 19. The multiplicity (by Cg) handles the gaps caused by points breaking the 
line alike condition with respect to the Euclidean metric. The role of the additive part, together with 
the second stage of the transformation, elaborated next, is to overcome the absence of points outside the 
path. In the second stage, for every pj with 1 < j < m, let 1 < j~ < j be the minimal index for which 
there exists u € T{pj-) with \pju\ < Ck • w{T(pj)), and let j < < m be the maximal index for which 
there exists u G T(pj+) with \pju\ < Ck ■ w(T[pj)), direct the sub-path between p^- and Pj+ towards pj. 
See Fig. for illustration. 



Figure 5: An illustration of the second stage in the adjustment transformation g. 

The indexing of points in Pm and notations introduced in the algorithm are used throughout this 
section. The tree T{pi) G i?, for 1 < i < z, is sometimes denoted by Ti for short. 

3.3 The Validity of the Output 

We consider each solution separately and show it forms a valid assignment p. 

Validity of Solution (i). Follows from the validity of the Hub algorithm 

Validity of Solution (ii). The subgraph of Gp induced by the points of Pm is strongly connected due 
to the validity of the Hub algorithm. All trees in R are bi-directed trees sharing a common point with 
Pm and therefore, the whole graph, Gp, is strongly connected. 

Validity of Solution (iii). Each tree in R induces a strongly connected subgraph of Gp, thus, it 
is suffices to show the connectivity of the minor obtained from Gp by contracting all trees of R. The 
subgraph of Gp induced by the points of the middle sub-path forms either one directed cycle or two 
directed cycles sharing common vertices. By the correctness of the Hub algorithm, each of the sub-paths 
Pi and Pr induces a strongly connected subgraph of Gp. In addition, each of them shares a common 
vertex with the middle sub-path, thus, the whole graph, Gp, is strongly connected. 

Validity of Solution (iv). Since we already verified the validity of Solution (iii), we are only left to 
show that each of P/ and Pr induce strongly connected subgraphs in Gp. We consider the sub-path Pi, 
while the case of Pr is symmetric. 

Let Pi : Pi ^ R+ denote the assignments obtained by applying the ID algorithm on the sub-path 
Pi with respect to hs. Due to the validity of pi, the graph induced by pi with respect to hs is strongly 
connected. Let {pi,pj) be an edge in this graph, we show that there exists a directed path from pi to pj 
in Gp. Assume w.l.o.g., 1 < i < j < 1. By the definition of hs, there exist u G T{pi') and v G Tipy) with 
1 < i' < i < j < j' < m, such that pi(pi) > hs(pi,Pj) = |nn|. Since the final assignment p is obtained 
after applying the adjustment transformation g, we have p{pi) > Cs|un|. 

Case 1: \uv\ < {ck - 1) • w{T{pii)) or \uv\ < {ck - 1) • w{T{pji)). 

Assume, w.l.o.g., that the second condition holds, then \pj'u\ < Ck ■ w{T{pj>)) and thus, by the 
definition of transformation g, the directed path from pu to pj’ and the edge {pj',u) are contained 
in Gp. Together with the directed path in T^pu) from u to the root pu they form a cycle containing 
both Pi and pj . 
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Case 2: both |uii| > (cfe — 1) • w{T{pii)) and |ui>| > (ck — 1) • w{T[pji)). 

Since G Pv, and Pi' is the output sub-path after performing the flatten procedure, then 

^Pi'iPiGPj') ^ Cs\pi'Pj'\. Therefore, we have 

\PiPj\ < 6p^,{pi,pj) = 6p^,{pi',pj') - {6p^,{pi',p,) + Sp^,{pj',pj)) 

< Cs\pppj'\ - {\piPj\ - {w{Ti') + w{Tj,) + |mi;|)) 

< cflwflTp) + w{Tj>) + \uv\) - \piPj\ + w{Ti') + w{Tp) + \uv\ 

< {cs + l){w{Ti') +w{Tj')) + {cs + l)|Mr'| - \piPj\ 

< {Cs + l){2\uv\/{ck - 1)) -I- (Cs -I- 1)|UW| - \piPj\ 

< (cs + 1)(2 |m?;|/(8(1 -|- Cs))) -|- (cs -|- 1)|m?;| - \p^pfl 

=> \PtPj\ < = Cs\uv\. 


3.4 The Approximation Ratio 

Let SOL denote the cost of the output of the algorithm for the input set S and let p* : S ^ ]R,+ denote 
an optimal assignment for S of cost OPT. We show that SOL < (1.5 — e)OPT. 

Let W = w{MST{S)), r = w{R)/W, cm denote the longest edge in MST{S) and I = 'w{eM)lW. 
As shown in [5], OPT > W + w{eM) = W{1 + 1). Next we show several upper bounds on the ratio 
SOL/OPT, corresponding to the four solutions computed during the algorithms and hnally conclude 
that the minimum among them equals at most (1.5 — e). 

Approximation bound for Solution (i). Due to the analysis of the Hub algorithm done in [5], we 
have SOL <W + {w{Pm) + w(eM))/2 = 1T(1.5 — r/2 + 1/2). Therefore, 

SOL ^ Wil.5-r/2 + l/2) 1.5 - (r - Z)/2 

OPT - W{1 + 1) “ 1+1 ■ 


Assume SOL > (1.5 — e)OPT, then 


1.5-e < 


SOL 

OPT 


< 


1.5- (r-Z)/2 

m 


< 1.5-r/2. 


implies r < 2e and the following corollary follows. 

Corollary 8. One of the following holds, SOL < (1.5 — e)OPT or r < 2e. 

The following lemma is crucial for introducing the bounds for Solutions (ii) and (iv). 

Lemma 9. Let Cg and the notation Qp be defined as in procedure f. Given a constant 5, 

1 . there exist two pairs of points {u,v), (y,w) connecting two disjoint sub-paths in Pm, each pair with 
stretch factor greater than Cg, that satisfy 5p^{u,v) > 6 and 5pj^[y,w) > S; or 

2 . there exists an edge e G Pm defining Pp and P^r (the two paths of Pm\&, to the left and to the 
right of e, respectively), such that for every {u,v) G Qp i UQp^r, 5pm{u,v) < S. 

Proof. If the first condition holds, we are done. Otherwise, fix e to be the rightmost edge in Pm- If the 
second condition does not hold for e, then there exists exactly one pair (u,v) G Qp i with (5p„ (tt, v) > 5. 
Replace e with the consecutive edge to its left in Pm- Continue the process until for every {u,v) G Qp i, 
Sp^j{u,v) < 6. Note that this condition hold when Pp contains a single edge. If the process ends with 
a separating edge e for which there exists a pair {u,v) G Qp,,r with 6 p^{u,v) > 5, then u is a common 
endpoint with the preceding edge in the process and there exists a point w G Pp such that the two pairs 
{w,u), (u,v) satisfy the first condition. □ 

Approximation bound for Solution (ii). Consider the value ChW, where Ch = 20e. One of the two 
conditions of Lemma denoted by 1^1 and L|^ 2 , respectively, must hold for <5 = ChW. We start by 
assuming that condition ij^l holds which leads to Lemma [lOl 

For every q G Pm, we have that the cost of Solution (ii) equals at most w{Pm) + niax{|gpi|, |gpz|} -I- 
2w{R). Consider the path ^f{PM) obtained from Pm after applying the flatten procedure / on the sub¬ 
paths Pgi and Pe'-. Note that the length of this path is at most w{Pm) — 2chW{l — l/cg) and it shares 
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common endpoints with Pm- Let c be the point on ^/{Pm) closest to its midpoint. The midpoint may 
lie on an edge of Pm or on a shortcut performed by /. If it lies on a shortcut, we undo it and only one 
shortcut remains. Thus, the point c is at Euclidean distance at most ^[w{Pm) — ChW{l — + 'w{eM)], 

from both endpoints and we have 

SOL <W{l-r+l[{l-r)- c„(l --) + 1] + 2r) = W[1.5 -- I - r)]. 

This implies, 

SOL ^ 1-5 - - I - r) 

OPT - 1 + 1 ■ ^ ’ 

Lemma 10. If condition holds for 6 = ChW then SOL < (1.5 — e)OPT. 

Proof. Assume towards contradiction that SOL > (1.5 — e)OPT, then by Corollary]^ r < 2e and to¬ 
gether with equation ([^ we receive 


1.5 - i(c/j(l - - Z - r) 1 1 

1.5-e< -——^- < 1.5--(cft l--r) 

1 -|- 1 2 Cs 


< 1.5- J(20e(i)-2e) < 1.5 - e 

Z 0 


contradiction. 


□ 


From now on, assume condition 1^2 holds for S = ChW and let e G Pm be the edge satisfying the 
condition. Let t = {J2(u,v)GQr + I2(u,v)eQp , / W, we give an additional bound 

to the cost of Solution (ii) using the same arguments used for the case where condition IMl holds. 
Let c be the point on "^/(Pm) closest to its midpoint. Since every pair (u,v) G Qp^i U Qp^r satisfies 
5p„{u, v) < ChW, the point c is at Euclidean distance at most ^[w{Pm) — {t — Ch)-W{l— -l-■lc(eM)] 
from both endpoints, thus 

S'OL ^ 1-5- 5 ((t-c/,)(l- 

OPT - 1 + 1 ■ ^ ’ 

Note that although the analysis for equation considers the path ^f{PM), the flatten procedure / is 
performed after Solution (ii) is computed. 

In the analysis of Solution (hi) and (iv) we have w{R) < (r -|- t)W , since applying / on Pgi and 
Pgr- moves portions of the paths connecting points in Qp , U Qp^^ to R. Consider the iteration of the 
algorithm for the edge e and a choice of 4 points pi , pi> , pp , Pr G Pm satisfying: pr is the rightmost point 
in Pm with Ur G T(pfl) connected (at any direction) in Gp* to a point in T{pi>) for pr to the left of 
e and, symmetrically, pi is the leftmost point in Pm with ui G T{pi) connected (at any direction) to a 
point in Tfpp ) for pp to the right of e. Meaning, there is no edge in Gp* connecting between a point in 
T{p) for p G P/\{p/} and a point in T{q) for q G Pp and no edge connecting between a point in T{p) for 
p G Pr\{Pr} and a point in T{q) for q G Pv . Let x denote the ratio w{Px)/W. 


Approximation bound for Solution (iii). Preforming the Hub algorithm on Pi and Pr, separately, 
result in two assignments with a total cost of at most 1.5(1T — w{Px) — w{R)) + w{eM)/2 + w{eM)/2. 
Directing the path Px and the trees in R, and assigning all roots in R their assignment, together 
with adding the two edges, {ui,up) for ui G T[pi),up G T{pp) of minimal length and [uipUr) for 
Ur G T{pii),Ur G T{pr) of minimal length (or the edge {pi,Pr) instead if it is cheaper) costs at most 
w{Px) + 2 • w{R) + \uiup\ + \uiiUr\. Overall, we have a total cost of at most W + ^(1 — a: -|- (r -|-1)) • 
W + IW + \uiup I -I- \urUr\. 

Since there is an edge connecting a point in T{pi) with a point in T{pp) (in some direction) and 
an edge connecting a point in T{pit) with a point in T{pr) in the optimal solution, we have OPT > 
W + \uiUr' I -I- \ui>Ur\ — w^cm), hence. 


SOL IT (1 — / -t-^(l — X + {r -t-t))-t-21)-t- \uiUr' \ + \uiiUr I 

OPT ~ W{1 — 1) + \uiUr'\ + \uiiUr\ 


< 1 + 


W{^{1 - X + {r + t)) + 21) 
W{1-1) + \uiUp \ + \ui'Ur\ 


< 1+^(1 


a: “t" (r -t- f)) + 21. 


(4) 
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Approximation bound for Solution (iv). Let pi : Pi ^ R'*' and pr '■ Pr ^ R’*' denote the 
assignments obtained by applying the ID RA algorithm on Pi and Pr, respectively, with respect to the 
distance function hs- Let p' : Pi U Pr ^ R^ denote the union of the two assignments and let OPT' 
denote the cost of p', i.e., OPT' = p'(''’)■ 

Claim 1. OPT' < OPT. 

Proof. We show that the optimal assignment p* can be adjusted to an assignment 
p : P/ U Pf —t R"^, valid for Pi and P^, separately, with respect to hs, of the same cost. 

We define, 

p{pi)= max {p*(u)}, p(pr-) = max {p*(u)}, 

vGT{pi), vGT{pi), 

l<i<m 7n-\-l<i<r 

and for every pj with j G {1, l-l}U{r + l,z}, p{pj) = maXr^Tipj){p*{v)}- 

Let u and v be two points in the same sub-path, w.l.o.g., Pi, and let {u = ui,U 2 , .■.,Uk = v) be the 
path from it to u in Gp*. Consider the sequence (u = yi,y 2 , ...,yk = v), obtained by replacing every 
Ui € T{Pi) with i/i = r(ui), and every Ui G T{P„; U Pr) with p, = pi, for 1 < i < fc. We prove the above 
sequence forms a path from it to u in the graph induced by p with respect to P/ and hs, and conclude 
that p is valid and OPT' < OPT. Consider a pair of consecutive nodes in the above sequence, (yi, yt+i). 
Note that if Ui = pj (resp. itt+i = pj) for 1 < j < I, than Ui+i = ph (resp. Ui = ph) for 1 < h < m; thus, 
for every 1 < i < k, hy the definition of p and hs we have, p{yi) > p*{ui) > |itiiti+i| > hs{yi,yi+i) and 
the claim follows. □ 

Let ^g{p') '. Pi'UPr' -G R+ be the union of the assignments obtained after applying the transformation 
g on each of P; and Pr, separately. The following lemma bounds its cost. 

Lemma 11. cost[^g{p')) < [Cg -I- (c^ -I- 2cs(cfc -I- l))(r -|- t)]OPT. 

Proof. Consider applying transformation g on Pi and Pr- By multiplying the range of every point in 
Pi U Pr by Cg we obtain an assignment of cost Cg • OPT'. As for the additive part and the second 
stage, we analyze the cost with respect to Pp, while the case for Pr' is symmetric. Every pj G Pr, is 
responsible for an additional of at most Xj = Ck ■ w{Tj) + ,Pj+) < Ck • w{Tj) + Cs\pj-Pj+\ = 

Ck • w(Tj) -I- Cg[2cfc • w{Tj) + w{Tj-) + w{Tj+)] to the total cost, where Pj- and Pj+ are defined as in the 
definition of g. The first element in the summation is the range added to pj itself, and the second is the 
cost of directing the path between pj- and Pj+ towards Pj (depicted in Fig. [^in black). 

Allegedly, for computing cost{^g{p')) we should sum Xj over all pj G Pr U Pp and add it to the cost 
Cg • OPT', however, we observe that it is suffices to consider a point pi only once as Pj -, for the rightmost 
point Pj such that i = j~ and only once as Pj+, for the leftmost point pk such that i = k'^. Thus, we 
can charge Pj- and Pj+ themselves once on each of the elements Cg • w{Tj~) and Cg • w{Tj+) in the overall 
summation. Namely, charge every point pj for a total range increase of Yj = w{Tj)[ck + Cg(2cfc -I- 2)]. 
Summing Yj over all pj G Pr U Pp and adding the cost Cg • OPT', using Claim gives 

cost{g{rho')) < Cg • OPT' + w{Tj)[ck -f Cg(2cfc -I- 2)] 

= Cg • OPT' + [ck + 2cs(cfe -I- l)]{w{R)) 

^ [cs + (cfc -l- 2cg(cfe -|- !))(?’ -f t)]OPT. 


□ 


Note that g has already assigned to every pj G Pr U Pp an assignment greater than w{T{pj)). 
Directing all trees in R towards their roots, directing the path P^ and adding the edge {pi,Pr), adds to 
the cost at most (2x -I- (r -|- f))W < {2x -I- (r -|- t))OPT, and together with Lemma 11 we receive 


SOL 

OPT 


^ Cg -|- {ck 2cg(cfe -|- 1) -|- l)(r -|- t) -|- 2x 


(5) 


Lemma 12. If condition ^^2 holds for d = ChW, then SOL < (1.5 — e)OPT. 
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Proof. Assume towards contradiction that SOL > (1.5 — e)OPT. By Corollary r < 2e and together 
with equation ([^ we receive, 


1.5 — i((t — Cft,)(l —— I — r) 1 1 

1.5 - e < - il, ^ < 1-5 - 2 - 20^)5 - 


<1.5-A+3e 


t < 40e. 


Replacing r and t with the above upper bounds in equation ([^ gives, 

1.5 — e < Cg + (c/g + 2.Cg(^Cf^ + 1) + l)(r + t) + 2x 

< 1^+70(42e)+ 2x ^ x > ^ - 1472e, 

4 ^ ^ 8 

and by equation Q we have, 

1.5 — e<lH-—(1 — X {r 1)) + 2/ < 1 + —(1 — (— — 1472e) + 42e) + 21 

2 2 o 

< 1.5-4+ 757e + 21 ^ / > 4 - 380e. 


The upper bound on I together with equation 0 imply, 

1 


1 

1 + 1 2 


1 

1 + Z ^ 2 


1 + ^ - 380e 


e > 


8 

105 


in contradiction to our choice of e = 


5 

10 ^ • 


We conclude with Theorem [T^ derived from Lemma together with Lemmas 10 and[l^ 


□ 


Theorem 13. Given a set S of points in for d > 2, a minimum cost range assignment (1.5 — e)- 
approximation can be computed in polynomial time for S, where e = . 


The reader can notice that our algorithm yields a better approximation bound than stated in the 
above theorem. However, we preferred the simplicity of presentation over a more complicated analysis 
resulting in a tighter bound. 
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